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Abstract. We build a simple model for feedback systems involving small RNA 
(sRNA) molecules based on the iron metabolism system in the bacterium E. coli, and 
, compare it with the corresponding system in H. pylori which uses purely transcriptional 

^H' regulation. This reveals several unique features of sRNA based regulation that could 

be exploited by cells. Firstly, we show that sRNA regulation can maintain a smaller 
I turnover of target mRNAs than transcriptional regulation, without sacrificing the speed 

^ ' of response to external shocks. Secondly, we propose that a single sRNA can prioritize 

the usage of different target mRNAs. This suggests that sRNA regulation would 
be more common in more complex systems which need to co-regulate many mRNAs 
efficiently. 
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! 1. Introduction 

Small noncoding RNAs (sRNAs) have recently been discovered as key components of 
genetic regulation in systems ranging from bacteria to mammals [H [21 [3], and this 
has spurred much activity in understanding their functional advantages. For example, 
degradation of sRNA with their target mRNAs has been proposed as a mechanism to 
obtain ultrasensitivity [1]. 

One particular system which has been extensively studied is the Fe-Fur system in 
the bacterium Escherichia coh which contains the regulatory sRNA RyhB [5], [6]. This 
system is responsible for maintaining homeostasis of Fe"^"*" ions, which are essential for 
cell functioning but also poisonous at high concentrations. During aerobic exponential 
growth, iron-using enzymes in E. coh utilize around 10^ Fe atoms per cell generation, 
but more than 10"^ Fe"'"+ ions in free or loosely-bound form is poisonous [3 [8]. Thus, the 
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Figure 1. (a) The sRNA motif found in E. coli, consists of three variables, /, the 
Fe-Fur complex which depends on the amount of loosely bound iron, r, the sRNA 
ryhB, and m, the mRNA of iron-using proteins. The red barred arrows between r and 
m represent the formation of the r — m complex and its subsequent degradation. The 
other red barred arrow indicates transcriptional repression of r by /. The influx of /, 
denoted by C, is divided into the channel A, regulated by m (green arrow), and the 
non-regulated channel B. (b) The transcription motif found in H. pylori. Here, m is 
transcriptionally activated by /. 

cell faces the problem of maintaining a huge flux of Fe through a small reservoir, and 
at the same time channeling this flux into its most essential functions when the cell is 
starved of iron: It needs to prioritize and sort the usage of a limited resource. 

Reference [9] describes a detailed model of the Fe regulation in E. coli, which 
incorporates several feedback mechanisms that together secure the system against both 
up and down shifts of the iron level. In this paper, we focus on the prioritization of the 
usage of iron by the sRNA regulation and its role in sudden iron depletion. The model 
in [9j is simplified into a core "motif (Figure [TK) describing the negative feedback 
used in iron homeostasis. The motif consists of three variables, / the iron-activated Fur 
(Ferric uptake regulator) protein complex that senses the Fe"*"^ level, r, the sRNA RyhB, 
and m, the mRNA of iron-using proteins. The sRNA works by binding strongly to the 
mRNA, after which this entire complex is rapidly degraded [TOt [111 [12]. The sRNA is in 
turn transcriptionally repressed by /. Thus, / effectively activates m, through a double 
negative link via the sRNA. 

Interestingly, the regulation of iron homeostasis in the bacterium Helicobacter pylori 
differs from that of E. coli in one important respect: the regulation via the sRNA 
RyhB is replaced by a direct transcriptional activation of the mRNA m by / (see 
Figure [lb) [I3l [H] . These two bacteria motivate us to compare the "sRNA motif with 
the corresponding "transcription motif . 

By studying these motifs, we demonstrate the following interesting aspects of sRNA 
regulation: (i) both motifs can be adjusted to have similar response times, but the 
metabolic cost is different for the two motifs, (ii) a single type of sRNA can efficiently 
prioritize expression level of various downstream target mRNAs, thus prioritizing the 
usage of a limiting resource. We also discuss possible experiments to test these results. 
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2. The Two Motifs 

2.1. sRNA moUf 

In the motif of Figure [1^ we focus on the case where sRNA binds to mRNA and both 
RNAs in the complex are degraded. In terms of an effective rate constant for the overall 
degradation, S, the dynamics of the concentrations of the sRNA (r) and its target mRNA 
(m) can be described by 

dr r 

— = OLr ■ r ■ m (1) 

dt Tr ^ ^ 

dm m , , 

— = am d-r -m (2) 

Tmrna 

where and set the respective production rates, and and define the 
background degradation times. We next rescale the parameters by measuring 
concentrations in units of am^r and measuring time in units of Tj. 0. 

To these rescaled equations we also add an equation for /, a small-molecule- 
activated transcription factor for the sRNA (as shown in Figure [1^). We simplify this 
two step reaction from / to regulation by assuming that all bindings are first order, and 
by rescaling binding constants such that / = 1 results in half-repression of the promoter 
of the sRNA gene. Including the import and consumption of / in analogy to the Fe-Fur 
system[9], the full dynamics of the motif becomes: 

# ^ I C-A-m -B-f when / > 0, 

dt \C when / = 0, ^ ^ 

dr a , ^ , 

-=^-r-7r-m, (4) 

dm m 

— = 1 - — - 7r ■ m. (5) 

Lib Tffi 

Here the dimensionless mRNA degradation time = TmmalTr-, and the dimensionless 
degradation rate of the RNAs is related to the dimensionfuU parameters as follows: 

7 = bararl- (6) 

In the absence of m, r is degraded relatively slowly because unbound sRNA are 
quite stable in vivo [llj. This means that r,., which is our rescaled time unit in eqs. 
dS])-©, is around (25/ In 2) min. In the absence of r, m is degraded with = 0.2, the 
lifetime of RyhB target mRNA in E. co\\ from ref. pLlJ. 

The most important parameters that determine the dynamics of sRNA regulation 
are a and 7 in eqs.Q and ([5]). We can determine a reasonable range of values for 
these two parameters using experimental data. First, for the Fe-Fur system a and 7 
are mutually constrained by the observation that the target mRNA, sodB, is depleted 

X The equation are formally rescaled by replacing m — > m/{amTr), r — > r/{amTr) and t — > t/r^. 
Thereby, a unit of time corresponds to the degradation time of the sRNA, and the unit of production 
rates corresponds to the production rate of the target mRNA. 
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Figure 2. Limits on rescaled sRNA-target mRNA coupled degradation rate 7 = 
6amT^- The colour shading indicate that the model takes longer (shorter) time in a 
darker (yellow) region to be depleted 5-fold after the small RNA is activated by setting 
/ from 40 to zero at t = 0. The green solid line shows the contour line of 3 min. When 
production of mRNA is much faster than production of sRNA the 5-fold reduction can 
be obtained by a small 7, whereas a relatively small a = am/ctr makes production of 
sRNA so slow that the mutual degradation must be very fast. The response time is 
calculated based on the degradation time of sRNA = (25/ln2)min {9]. 



around 5- fold within 3 minutes after full induction of the RyhB promoter pjj. These 3 
minutes include the time required for RyhB production, set by a, plus the time for the 
produced RyhB to bind to sodB message and to degrade the complex, set by 7. For 
a < 1 there will never be enough RyhB to deplete the message completely, whereas for 
a slightly higher than 1 an extremely high 7 is needed for efficient depletion. 

Fig. [2] illustrates this. We simulate eqs. (jlj) and (ISl) to measure how long it 
takes for m to be depleted 5-fold after fully-activating sRNA at t = by changing / 
manually from 40 to zero. The solid line in Fig. [2] shows values of 7 and a that give a 3 
minutes depletion-time for sodB mRNA. With a values in the range from 3 to 10, we see 
that physiological 7 values would need to be between 100 and 10000. In addition, the 
reduction of iron consumption occurs as soon as the sRNA-mRNA complex is formed, 
even before the mRNA is degraded, but the measurement of [11] does not distinguish 
RNA species in complexes from the free form. Thus, the effective 7 could be larger than 
the above estimate. The value of a can be estimated from other data in ref. [H] of 
RyhB and sodB time series after induction and repression of RyhB. From these data, 
we estimate that a is between 2 and 5. In the rest of the paper, we explore the small 
RNA motif for a range of a and 7 values, around the estimates made above. 

The three other parameters. A, B, and C, in the iron flux equation ([3]), we set 
using experimental data [9|, as described in the appendix. 
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Figure 3. (a) The steady state level of r/{a — 1) (thin lines) and m/Tm (thick solid 
lines) vs. / for a — A with 7 = 15 (red solid lines) and 7 = 1500 (green dashed 
lines) in the sRNA motif, (b) The steady state level of m/Tm (thick lines) vs. / for 
the corresponding transcription motif parameters (see text). The time evolution of 
/ (divided by 200, thin lines) and m (thick lines) for the systems in (a) and (b) are 
shown in (c) and (d), respectively, after the sudden drop of C at i = 0. 



2.2. Transcriptional motif 

To model the transcription motif (Figure Wp) we replace (jl]) and (I5l) by the single 
equation: 

— = D - — (7) 

dt ' f^ + K^ Tm' ^' 

The first term models the direct transcriptional activation of m production by /, where 

D sets the maximum production rate of m and Kt sets the binding constant between 

/ and the DNA. The "Hill coefficient" h sets the steepness of the response, and is 

related to the cooperativity in binding. We use the same values of the influx C and two 

constants A and i? in ([3]) as those used in the sRNA motif. We choose D and in 

the transcription motif to have the identical steady state values of the / and m to the 

corresponding sRNA motif (parameterized by a and 7) for high iron (/ = 40) and low 

iron (/ = 5). It is not obvious that this is possible for all values of a and 7. In fact, we 

found that it is not possible for h < 2, but when h = 3, we can set D and Kt such that 

the above conditions are fulfilled for a G [0, 20] and 7 G [0, 2000]. Therefore, henceforth 

we keep h = 3. 
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Figure 4. The contour plot of sRNA motif response time, At, at which / recovers to 
and subsequently remains within 95% of the final steady state value. The corresponding 
transcription motif systems show a faster response in the shaded region. 

3. Results 

3.1. Degradation times 

Figures [3^1 and b show how the steady state values of r and m depend on the steady 
state / level, in the sRNA motif for two different 7 values, and, respectively, in the 
corresponding transcription motif systems. For the sRNA motif, a larger 7 results 
in a much larger ratio between maximum and minimum m values, and a steeper 
drop between them (Figures [3^). This is expected because in the 7 — > cxd limit, 
m = max(0(l/7), rm[l — a/(l + /)]), and this steepness in the sRNA regulation is 
referred to as "ultrasensitivity" in For the transcription motif, on the other hand, 
the dependence of m on / is weaker and nearly unaffected by the corresponding change 
in D and Kt (Figure [3)d). This is because the steepness of the curve in the transcription 
motif is determined by the Hill coefficient h, which is set to be relatively high value 3 
but still not enough to give as sharp slope as the sRNA motif with 7 = 1500. 

The importance of 7 is also evident in the dynamical response of the motifs to 
a sudden depletion in the external iron source C (Figures [3t and d). For the sRNA 
motif. Figure [3t; shows that a larger 7 results in a faster drop in m level, and a quicker 
approach to the new steady state level. The same is true for the / level also. That 
is, a large 7 naturally ensures a faster removal of all excess m, while also allowing / 
and m to climb back to non-zero steady levels even after / drops almost to zero during 
the initial shock. In contrast, the transcription motif displays approximately the same 
timescale of mRNA drop for the two corresponding cases because a drop in m takes a 
time proportional to r^, independent of other parameters. 

We investigate this further by quantifying the response time. Looking at Figure [St, 
we see that the / level drops very sharply and then rises towards its final steady state 
value. During this rise, at some time. At, f reaches within 95% of the final steady 
state value. We use At as a measure of the response speed of the system. In Figure H] 
we show a contour plot of At for a range of values of a and 7. The corresponding 
transcription motif systems show a faster response than the sRNA motif in the shaded 
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region. The comparison indicates that sRNA based translational regulation produces 
a faster response than transcriptional regulation when a > 2 (and a is not too large 
compared to the investigated range of /) and 7 is larger than a critical level, around 
150 (the unshaded region in Figure Hj). 



3.2. Metabolic cost 

The relatively slow response of the transcription motif is, of course, because its response 
timescale is set by Tm, which we keep constant. We emphasize that it is possible to 
achieve a fast response in the transcription motif also by decreasing Tm- There are some 
costs associated with this alternative strategy though. Faster mRNA turnover due to 
lower Tm requires a higher production to maintain the same homeostatic levels of /. At 
high / levels, where mRNA is high, the sRNA motif secures a low degradation rate of 
mRNA, whereas the transcription motif produces the mRNA at its maximum rate. On 
the other hand, at low / the sRNA motif maintains a high rate of mRNA degradation, 
while the transcription motif saves resources by reducing mRNA production. Therefore, 
for a given response speed to sudden iron starvation, the sRNA motif is less costly if 
the bacterium usually lives in iron-rich conditions, whereas the transcription motif is 
preferable if the organism mostly lives in iron-poor conditions. 



3. 3. sRNA prioritize downstream protein levels 

3.3.1. Prioritization in the iron feedback motif The properties of sRNA regulation can 
be used in an interesting way when more than one kind of mRNA is under regulation 
[T2] . Different mRNA can have hugely different binding strengths to the regulating 
sRNA, and thereby very different effective degradation timescales. This is because the 
binding strength of the r-m complex is, to a first approximation, an exponential function 
of the number of matching base-pairs, which can vary by an order of magnitude across 
different mRNA: The free-energy gain per matching base-pair is around 1 to 2 kcal/mol, 
which can give the difference in statistical weight exp^AG/ksT) ^ 5 to 30. Because of 
this property, sRNA regulation could be used to prioritize the degradation of different 
mRNA. We illustrate this by adding a second mRNA to the sRNA motif: 

- r ■ (7imi + 72m2) (8) 

- 71^1, (9) 

- 72rm2. (10) 



dr 


a 








— r 


'dt 


1 + / 


dm I 




mi 


dt 






dm2 




m2 


dt 




■ (ED, 


replacini 


g m 



For /, 

different effective degradation rates, 71 and 72, resulting in different steady state levels 
for a given C value. Other parameters are set as in the case with a single mRNA. 
Figure O shows the behavior of this system. The steady state level in Figure shows 
that m2, the mRNA with larger 7, is suppressed more than mi on depletion of /. 
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Figure 5. Investigation of mRNA prioritization by the sRNA motif, (a) The steady 
state level of r/(Q!—l) (blue dotted line) , mi/{Tmami) (red solid line), and m2 / {TmOirm) 
(green dashed line) vs. / for a = 4, 71 = 150, and 72 — 1500. (b) The time evolution 
of / (divided by 400, blue dotted line), toi (red solid line), and m2 (green dashed hne) 
after the sudden drop of C at t = 0. 



The prioritization behavior is easily understood by considering the extreme case 
where both 7s are very large but 72 S> 71. Then, depending on the level of /, the steady 
state is such that either (i) both mRNAs are near-zero, (ii) only mi is non-zero, (iii) 
both mi and m2 are non-zero. 

Taking into account finite 7 values, the prioritization efficiency becomes, 
respectively. 



(iii) mi ^ r^a„, and ma ^ r„(a„2 - j^) for large /, i.e. ^ - 1 < /• 

where the 0{x) are some functions that are small and proportional to x for small x. 

Case (ii) is what we refer to as the "prioritized state" , where only one of the mRNA 
(the one with smaller 7) is present and the other's level drops to near-zero. Clearly, the 
larger the difference between the 7 values of the mRNAs, the better the prioritization; 
the order of magnitude difference in 7 is important to have the clear prioritization. In 
addition, a, the sRNA production rate, determines the range of / for which the sRNA 
is affected, and a larger a results in prioritization being effective for a wider range of /. 

In addition, the dynamics in Figure [5]d shows that, when C is suddenly dropped, 
the mRNA with a larger 7 value is rapidly depleted while the other mRNA stays at a 
higher level. That is, the sRNA is not only able to prioritize the mRNA steady state 
levels, but is also able to remove the "unwanted" mRNA m2 much quicker than the 
"wanted" mRNA mi. 




1< 



3.3.2. Prioritization of multiple mRNAs It is clear that the prioritization can be 
generalized to the case of more than two kinds of mRNAs regulated by a single type 
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Figure 6. Prioritization of mRNAs by (a) sRNA regulation and (b) transcriptional 
regulation. The steady state level of mRNAs normalized by their maximum value 
is plotted against (3. (a) The prioritization by sRNA regulation, for ami — ^rn2 = 
ama = cxmi = 0.25, and 71 = 10^,72 = 10^,73 = 10"*, 74 = 10^. The value of ji 
are chosen to be large and have 10-fold difference between different mRNAs, so that 
the separation becomes clear, (b) Transcriptional regulation, where the transcription 
factor has concentration (3 and acts as a repressor. The Hill coefficient h is 3. The 
binding constants are Ki ~ 1, K2 = 0.75, K3 = 0.5, and K4 = 0.25. The values of Ki 
are chosen to be Ki — '^i^ th.at the value of /3 at which the mRNA starts to 

be significantly degraded is the same for the two motifs. 



of sRNA. To illustrate this we ignore feedback through / and consider the following 
general system with n different types of mRNAs: 

^ = /? -r - ^7irmi, (11) 

dmi mi , . 

—rr = ami lirrrii for t = l,n. (12) 

at Tm 

Here, to focus on the prioritization, the production term of sRNA in ([8]) which 
contains feedback from iron concentration is replaced by a constant (3. We assume 
7i < 72 < ■ ■ ■ < 7i < 7i+i < ■ ■ ■ < In without loss of generality, and rescale all variables 
to be dimensionless such that X^iLi '^rni and the degradation time of sRNA are unity. 

Figure shows the normalized steady state level of mRNAs, mi/{am,Tm), versus 
the production rate of the sRNA, (3. As (3 becomes larger, sRNA increases and more 
mRNAs are degraded. As shown in Figure [6k, the n-th mRNA with the largest 7 
is degraded first. This also "protects" other mRNAs from degradation because the 
sRNAs are also degraded together with the n-th mRNAs. The {n — l)-th mRNA 
starts to be degraded when the level of n-th mRNA becomes low enough, which 
occurs roughly at /3 ~ a„. The {n — 2)-th mRNA starts to be degraded when 
/5 ~ a„ + ctn-i, and so on, and finally all the mRNAs are almost completely degraded 
when (3 ~ ^"=1 = 1- The separation of the level between the [k + l)-th mRNA and 
the k-th mRNA foT_y2i=k+i < (3 < X^ILfc ^rni becomes clearer for larger difference 
between 7j+i and 7J§|. This multistep-switch-like degradation upon changing the value 

§ In the case 7^+1 /7i <C 1 for any i, the steady state level of mRNA for Yll^k+i '^m.i < P < X)"=fc '^mi 




Figure 7. Schematic picture of the prioritization of various kinds of mRNAs by a 
single type of sRNA. Different kinds of mRNAs are represented by different colors, 
and have different value of 7 (left boxes), (a) When a small amount of sRNAs is 
produced, the mRNAs with larger 7 are degraded (shaded region in the right box), 
while the levels of mRNAs with small 7 are scarcely affected, (b) As more sRNA is 
produced, the mRNAs with smaller 7 are also degraded. 



of (3 is the characteristic feature of the sRNA regulation. This prioritization mechanism 
is schematically presented in Figure [71 where the mRNAs are degraded in descending 
order of the value of 7. 

We note that it is possible to have different steady state of mRNAs in the 
transcriptional motif also. Suppose r is the concentration of the repressor with a 
production rate jS. Then the simplest model for the mRNA level is 

I -P-r, (13) 
= — -r for t = l,n. (14) 

with a Hill coefficient h and binding constants Kf, we assume Ki > K2 > ■ ■ ■ > 
without loss of generahty. Here, the variables are rescaled to be dimensionless such that 
Ki and the degradation time of r are unity. The normalized steady state level is given 
by rrii/ {amiTm) = T+jj^/Ky^'- The levels decrease as P increases with the slope determined 
by the Hill coefficient h, and the characteristic value of /3 where the rrii level becomes 
half of its maximum is given by Ki. Figure [Hb shows the separation of the steady state 
by transcription regulation with a relatively high value of the Hill coefficient h = 3 and 
various values of Ki. We see that the separation of the rrii level is not as sharp as with 
sRNA regulation, and rrii does not change much upon changing jS especially for large 

is estimated as follows: mi/(r,„ami) w 1 for i < /c, mk/{T.marnk) ~ iYli=k^rni - P)/amk, and 
mi/{Tmai) ^ 1 for 2 > fc + 1. 
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Ki. The sRNA regulation is more effective in the sense that the prioritization of the 
mRNAs is sensitive to small changes in (3. 

4. Discussion 

Our analysis pinpoints three features that are particular to sRNA regulation in a 
feedback system. First, as sRNAs act through degradation, the regulation can, in 
principle, be very fast, generating a near instant response. Second, as the sRNA motif 
uses a double negative link, instead of a direct activation regulation, it has a higher 
metabolic cost for conditions where downstream targets are repressed. Third, and most 
interestingly, it has the capability of efficiently prioritizing the usage of downstream 
target genes. 

It is essential for the prioritized expression of downstream targets that the sRNA 
as well as the mRNAs are degraded together after they form a complex. In the large 
7 limit, the effective prioritization occurs because the degradation of m2 interferes with 
the degradation of nii by sequestering r, leaving less unbound r to bind mi. Indeed, if 
the sRNA simply catalyzed the degradation of the mRNAs without itself being degraded 
(i.e., no 7jmj terms in ^ or ( fTTi) ). the "protection" of mi is also lost. The degradation 
of different mRNAs doesn't interfere as in the previous case, and thus the prioritization 
is less effective. 

The switch-like behavior due to the "ultrasensitivity" of the sRNA regulation [1] 
together with the prioritization suggested in this work opens the possibility of more 
sophisticated regulation of gene expression. In particular, if each mRNA is targeted by 
several kinds of sRNAs, the combinatorics allows one to realize various logic gates. For 
example, if different kinds of sRNAs can bind to the same part of the targeted mRNA 
to trigger the regulation, it behaves as an "OR" gate. Such an example is known in V. 
cholerae, where four different sRNAs regulate HapR and any one of them is enough for 
regulation^. It is also, in principle, possible that one mRNA has multiple binding sites 
for different sRNAs (JJ]. In this case, if binding of all the sRNAs is necessary to trigger 
the regulation, it realizes an "AND" gate. The concentration dependent prioritization 
can add more complicated functions to the logic gates. It is an interesting future issue 
to explore the possibility of combinatoric regulation by sRNAs. 

We have tested numerically that the results of this paper do not strongly depend 
on the specific form of ([3]) for /. As long as in- and out-fiuxes are large enough to allow 
a much faster response of / than of m, then m is the rate-limiting factor. Another thing 
to note is that we assumed that the Hill coefficient for the repression of sRNA by / is 
1, though we introduced a Hill coefficient /i = 3 in the transcription motif to achieve 
the sharp contrast in m at high and low /. If we put a Hill coefficient H > 1 in the 
sRNA motif by replacing the production term of r with a/{l + f^), the m vs. / curve 
becomes even steeper, which makes the response sharper. 

II The authors do not know of an estabhshed example, but there exists an mRNA that has muhiple 
binding sites for the same micro RNA in eukaryotic cells [lH [16] 
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5. Experimental tests 

The suggested prioritization possibility for downstream mRNAs invites experimental 
tests. One key quantity of interest is the degradation parameter 7, which in fact sets 
the efficiency of the whole sRNA regulatory system. For large 7, the sRNA regulation 
works as a step function: when production of sRNA is larger than production of 
downstream targets, these are instantly removed. Therefore it is essential to measure 
degradation times of downstream targets under various expression levels of the sRNA. 
These degradation times are tightly coupled to the steady state level of downstream 
mRNA, and could therefore be obtained from bulk measurements, using for example 
microarrays for both RyhB and downstream targets. For a given r (=sRNA) in the 
steady state, the downstream mRNA level is mj(r) = Q;^./(l/r^ + jir), and therefore 
the slope of mi(0)/mi(r) — 1 versus r gives ■jiTm- 

An experimental test of the prioritization capability of the RyhB system is to 
consider homeostasis, and compare wild type with any mutant where the RyhB binding 
part of a downstream target gene is highly over-expressed. The over-expressed genes 
should be constructed such that they do not produce proteins that can bind Fe, and 
thereby indirectly influence the frcc/loosely bound Fe pool. For a given low level of 
external iron, the prioritized usage implies that for mutants where the over-expressed 
gene has small 7 the expression levels of other genes would be almost the same as 
in wild type. For the remaining, there will be large changes associated to RyhB being 
depleted by the over-expressed gene with large 7. If the difference in 7 between different 
target mRNA is of an order of magnitude, i.e., enough to have clear prioritization, the 
influence of the over-expressed gene should be quite sharp. Further, our prioritization 
scheme implies that as the external Fe is depleted further, the number of downstream 
genes which influence homeostasis should diminish monotonically. 

6. Conclusion 

Using the well characterized homeostatic response system for iron in bacteria, we 
analyzed the pros and cons of sRNA versus transcription regulation. The investigated 
negative feedback motif of Fe — > FeFur proteins — > Fe brings out the functional 
similarities and differences between the two alternative strategies of regulating 
downstream targets. For sufficiently high Hill coefficients, transcriptional regulation 
can reproduce the same steady state behavior as the sRNA regulation. Further, both 
regulations can in principle provide a fast response to sudden decreases in externally 
available iron. However, their functional capabilities differ in two important aspects. 

• First, adjusting parameters to obtain similar response times, the sRNA motif 
results in more turnover of target mRNAs in iron-poor conditions, whereas the 
transcription motif results in more turnover in iron-rich conditions. 

• Second, the sRNA allows a prioritization of expression level of downstream targets, 
thus efficiently regulating the usage of a limiting iron resource. At the same time. 
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unwanted mRNA is degraded more rapidly. This observation fits with the fact that 
the transcription motif is found in H. pylori, a bacterium with a small genome and 
limited capacity for genome regulation, while E. coli, which has a larger genome 
and can benefit from fine tuning of mRNA levels, has an sRNA motif. 

Our analysis suggests new ways to analyze other systems where multiple sRNAs 
regulate a more complicated response, including for example quorum sensing [1]. There, 
mutual binding inhibition between the sRNAs and a central regulator (LuxR mRNA in 
V. harveyi [1], HapR mRNA in V. cholerae [17], and the translational regulator RsmA 
in Pseudomonas [H]), may allow signals from adhesion and host factors to differentiate 
the sRNAs and thereby their downstream targets. Overall, our analysis suggests that 
sRNAs or micro RNAs may allow a near-digital reorganization of cellular composition, 
an observation which concurs with their ubiquity in regulatory processes associated with 
development. 
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Appendix: Parameters in the flux equation 

We set the parameters A, B and C in eq. (3) using data on the uptake and usage of Fe 
in E. coli and if. pylori. These systems are characterized by a huge flux of free Fe^^ , 
where the fluxes C, A and B are so large that the pool is replenished about 100 times 
per cell generation. 

In more detail, the constant incoming flux, C, is partitioned into two channels, A 
and B, that are motivated by the two separate ways of using iron in the Fe-Fur system. 
Flow through channel A is regulated, i.e., it can be reduced during iron-starvation, and 
is proportional to the mRNA level, but independent of / when there is any substantial 
amount of / in the system pj^ . In contrast, flow through channel B is proportional to 
/. In the case of the full Fe-Fur system [9] the regulation of this flux by other proteins 
becomes important when extracellular iron increases suddenly, but here we focus on 
the regulation by sRNA and therefore keep flow through channel B unregulated. That 
is, our motif is designed to respond to depletion of /. Our "minimal" motif cannot do 
without this B channel because removing it results in robust oscillations, for which there 
is no evidence in E. coli [5j. 

^ When we simulate eq.([3]), we used the form df /dt — C — A-m-[f /{f ^ Kcut)] — B ■ f with Kcut — 0.1 
to avoid the numerical difficulty due to the singularity at / = 0. This expression agrees with eq.([3]) in 
the limit of Kcut 0. 
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The parameters in are set as follows [S] : For conditions where extracellular iron 
is plentiful, we demand that the internal Fe"'"^ level is such that / = 40 |9] in steady 
state, while at the same time the net in- and out-fiux per cell generation is approximately 
100 times larger, representing the fast turnover and high usage of Fe [5l |9] . In addition, 
we assume that the iron flux is equally partitioned between the A and B channels, 
because we found it best fulfills the homeostatic requirements in our full model [9]. 
These conditions set the value of two parameters C = 4572 (i.e. C = 114 x / [9]) and 
B = 57. The value of A is completely determined by the steady state level of m, which 
depends on the parameters in the regulation part of the motif (i.e., a and 7 in the 
sRNA motif, or D and Kt in the transcription motif). For iron-starvation we demand 
that / = 5 [9], keeping the value of B and A constant. This condition is achieved by 
reducing C, reflecting the reduction in extracellular iron. The C change needed to get 
/ = 5 is dependent on the values of the regulation part of the motif. Note that, for a 
given regulatory motif, the steady state value of / is a unique function of the influx C. 
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